String structures in driven 3D complex plasma clusters 
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Abstract - The structure of driven three-dimensional complex plasma clusters was studied ex- 
perimentally. The clusters consisted of around 60 glass microspheres that were suspended in a 
plasma of rf discharge in argon. The particles were confined in a glass box with conductive yet 
transparent coating on its four side walls. This allowed manipulating the particle cluster by biasing 
the confining walls in a certain sequence and direct imaging of the cluster. In this work, a rotating 
electric field was used to drive the clusters. Depending on the field frequency, the clusters rotated 
(10 4 — 10 7 times slower than the rotating field) or remained stationary. The cluster structure was 
neither that of nested spherical shells nor simple chain structure. Strings of various lengths were 
found consisting of 2 to 5 particles, their spatial and temporal correlations were studied. The 
results are compared to recent simulations. 
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Introduction. — A complex (dusty) plasma consists 
of a weakly ionized gas and small solid particles 1 . The 
particle size can be from tens of nanometers to tens of 
microns. The particles are usually negatively charged due 
to the high thermal speed of plasma electrons. In ground- 
based experiments, the particles can be levitated against 
gravity by a strong electric field present in the plasma 
sheath (pre-sheath) around an electrode or wall [2p] . Hor- 
izontally, like-charged particles are held together by weak 
electric fields naturally present in the plasma or by exter- 
nal confinement. The latter can be provided, e.g. by a 
recession in the electrode [4] or by additional containers 
such as a glass box [5j|7] . 

The interparticle interaction in a suspended particle 
cloud consists of two parts: direct screened- Coulomb 
(Yukawa) repulsion and ion- wake mediated attraction. 
The ion wakes are formed by ions streaming past a micro- 
particle due to their focusing behind it pj-fll]. 

The structure of a suspended particle cloud is therefore 
determined by the interplay of the isotropic Yukawa-type 
and anisotropic wake-field forces giving rise to different 
competing symmetries and patterns, e.g. nested spheri- 
cal shells |5| and flow-aligned strings 12 13 . The order 



be governed by the ion flow Mach number in the dust dy- 
namics simulations [14] . 

Structuring is one of the basic processes in complex 
plasmas, as well as many other systems. For instance, 
dust chains, filaments, and stationary structures have been 
observed in sparks [l5] and dc plasmas [l6 17 . Note 



that competing repulsive-attractive interactions often de- 
fine the actual variety of observed strings, filaments and 
patterns (lattices) formed by the interacting strings as, for 
instance, in colloids 
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19 , in elec- 



in superconductors 
trons on superfluid helium [20], and in current- carrying 
plasmas [21] . 

To probe the structure of clusters in various states, some 
method of external manipulation is necessary. One way of 
manipulating a cluster is applying a torque to it. This can 
be done by means of external magnetic field [I], electric 
field 7,22 , or by the drag force of rotating neutral gas [9]. 

In the present paper, we studied a particular class of dy- 
namically driven clusters (DDC) - 3D clusters of charged 
particles driven by rotating electric field. The presence of 
the driving force deeply affected the actual cluster geome- 
try and the global symmetry of the cluster that appeared 
to be defined by spontaneously formed particle strings in- 
transition between these two configurations was found to teracting with each other. The recognition of string struc- 
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tures observed in DDC's is the main goal of this study. 

The lattices as well as competing symmetries shown by 
the string structures formed in our experiments resem- 
ble to some extent the Faraday wave patterns parametri- 
cally excited in water by vertical vibrations 23 . Similar- 
ities include, for instance, the competition between sub- 
harmonically generated equilateral triangles and regular 



neutral density filter camera B 



hexagons 24 and the super-lattice patterns [25]. There- 
fore, experiments with complex plasma clusters might help 
to get insight into generic structuring processes at the ki- 
netic level. 

Experimental setup. — Our experimental setup has 
been previously described in Ref. [22]. We used a ca- 
pacitively coupled rf discharge in argon at a pressure 
of 4 Pa. The inner dimensions of the chamber are 
14 cm x 14 cm x 26 cm [26]. The lower electrode diameter 
is 7 cm. The upper electrode is mounted at a distance of 
7.2 cm and has a diameter of 9.6 cm. The lower electrode 
is powered by a rf generator operating at 13.56 MHz and 
120 V peak-to-peak amplitude. The upper electrode and 
the surrounding walls are grounded. In the center of the 
upper electrode a grid-covered hole is made through which 
particles can be transferred into the chamber. The parti- 
cles used in these experiments are hollow glass spheres Q 
with a diameter of 22 ± 2 /am. 

On top of the lower rf electrode we placed a glass box 
as described in Refs. [7,22 . The inner surfaces of the box 
side plates are coated with Indium Tin Oxide (ITO) and 
are therefore conductive yet transparent. The plates are 
insulated from each other and from the lower electrode 
by plastic poles. Each plate is electrically connected to 
a separate function generator. The function generators 
output sinusoidal signals with a peak-to-peak amplitude 
of 20 V and frequency in the range of 100 Hz to 100 kHz. 
It is possible to make the resulting electric field inside the 
box to rotate in the horizontal plane by setting the phase 
shift of the sinusoidal voltage on the n-th plate to ±nir/2. 
This "rotating wall" technique was introduced in Ref. m 
to drive 2D clusters of particles. 

To create a cluster, microparticles were suspended in- 
side the glass box where they formed a spheroidal cloud 
with a diameter of ~ 5 mm. Applied rotating electric field 
contracted the cluster and forced it to rotate (in fact, co- 
rotate in the present experiments). The speed of rotation 
depended on the applied field frequency [7,22 . Cluster 
rotation was the fastest at 5 kHz and practically absent at 
1 kHz [22], therefore in this paper we used these two fre- 
quency settings to study rotating and stationary clusters. 

To image the particles, we used a 3D diagnostic method 
that is called stereoscopic digital in-line holography [26] . 
This method provides the information on 3D fully-resolved 
particle dynamics. A sketch of the laser and camera setup 
is shown in Fig. [I] The cameras operated at a frame rate 
of 50 frames per second during a time interval of 10 s. 
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Fig. 1: Sketch of the experimental setup. The microspheres are 
confined in a glass box placed on the lower electrode in argon rf 
discharge. 3D diagnostic system [26] consists of a laser (wave- 
length of 532 nm), beam splitter, and two identical channels 
consisting of a beam expander, mirror, neutral density filter, 
and digital camera. Interference patterns created by particles 
are registered directly by camera sensors. The resulting images 
contain complete 3D information on particle positions. 



No lenses were mounted on the cameras. The laser light 
was therefore registered directly by the camera sensors, 
after passing through the particle cloud. In the result- 
ing images, each particle was represented by a system of 
concentric interference rings, the depth information was 
encoded in the inter-ring spacing. 

Particle cloud parameters. — The hollow glass 
spheres used in our experiments are made of Scotchlite S22 
glass with a density of 2.5 g/cm 3 and have a wall thickness 
of 0.3 jam. They are sorted to achieve a narrow distribu- 
tion of particle diameters of 22 ± 2 /am. It is straightfor- 
ward to estimate the particle mass m = 1.11 ±0.23 ng 
and the neutral gas damping rat^7 e ff = 17.5 s _1 (at the 
argon pressure of 4 Pa). 

In our experimental conditions, the plasma parameters 
are estimated as follow^] the electron temperature is 
T e « 3 eV, electron-to-ion temperature ratio T e /Ti w 100, 



plasma density n ~ 10 cm 



The electron Debye screen- 
w 400 jam and the ion mean 



ing length is estimated as \]j e 
free path as X ia « 770 /am^\ 

The particle charge number is expected to be in the 
range of Z = 40 000 — 50 000, as standard theoretical ap- 
proximations allowed us to estimate (for instance, Z ~ 
46 000 follows from OML [29], Z ~ 39 000 from DML [30], 
or Z ~ 41000 from modified OML [3l]). Note that in 
our experimental conditions the ratio Np/Z ~ 2-3 is not 
large, where Njj is the average number of electrons inside 
the Debye sphere. Therefore, the actual particle charge is 



expected to be smaller due to electron depletion 32 

In experiments with clusters of only 2 — 3 charged 
particles (arranged as a string), the particle charge and 



1 "Glashohlkugeln" or "glass bubbles" produced by 3M GmbH 
Deutschland, IndustriestraBe, 55743 Idar-Oberstein. 



2 When calculating j e ^, one has to take into account that the 
particles are hollow, unlike discussed in Ref. [27| . 

3 Typical bulk plasma parameters, see Ref. [22] . 

4 Ion- atom charge exchange collisions [28] are assumed to be dom- 
inant. 
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damping rate can be measured directly by exploring the 
response of a cluster to harmonic excitation applied to 
the lower electrode [33] . A series of such measurements 
showed, for instance, that in doublets the upper particle 
charge number was about Z = 40 000, while the lower par- 
ticle's charge was up to 30% smaller. The effective damp- 
ing rate (~ 10 s _1 ) was measured to be in a reasonably 
good agreement with the above predictions. 

Given the approximate values of particle charge and 
mass and assuming that the simple balance equation mg = 
ZeE is valid, where g is the acceleration of gravity, e is 
the elementary charge, and E is the local electric field at 
the particle levitation height, one can estimate the Mach 
numbeij^] of the ion flow. In the mobility-limited regime 
the following relationship can be used: 



M 



mgX ia 
ZT P 



(1) 



It yields M = 0.6 — 0.7 in our experimental conditions. 
This finding is important for the analysis of cluster struc- 
ture in the next section, since it allows us to compare our 
experimental results to the dust dynamics simulation re- 



ported in Ref. 14 , where M was found to govern an order 



transition in 3D clusters. 

Competing symmetries in driven clusters. — To 

identify string structures in particle clusters, we analyzed 
the spatial and temporal correlations in particle posi- 
tion^ 

Spatial correlation. Horizontal distances between all 
particles in the (x,y) projection (top view) were cal- 
culated, resulting in a distance matrix. Particles were 
considered as coupled in a string if the horizontal dis- 
tance between them was below the threshold introduced 
by the "critical cylinder" (oriented vertically) of radius 
r cy i = 200 fim « ^A^ e , as illustrated schematically in 
Figs. [2ja),(b). (As it will become clear later, the chosen 
value of r cy i is smaller than the mean interstring distance 
in the (x, y) plane.) 

The top view of the cluster in Fig. [2|a) indicates that 
the majority of particles are organized in vertical strings. 
The side-view cross-section in Fig. [2^b) also clearly shows 
particle strings. Only few particles remain alone. Fur- 
thermore, the spacing between two strings is larger than 
the characteristic horizontal extent of a string. Therefore, 
most particles are well aligned in the vertical direction 
(along the z-axis). 

Temporal correlation. Another important criterion is 
the lifetime of a particle string. Studying the lifetime of 
particle associations one can easily distinguish pairs of par- 
ticles which only occasionally happened to be close to each 
other from persistent particle pairs or strings where parti- 
cles interact with each other by the wake-field attraction. 



3 Here, the ratio M 



Ui/c s 



of the ion flow velocity ui to the 
,-, where ra? is the ion mass. 



speed of ion sound c s i = y/Te/' 

6 The details of computational procedure will be published else- 
where. 
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Fig. 2: String structure of a rotating 3D cluster, (a) In the top 
view, the (x, y) particle positions (black dots) are structured ei- 
ther as a system of cylindrical shells (dashed lines) or tetragons 
(dotted lines). The inner dashed circle has a radius of 0.85 mm, 
the outer of 1.7 mm. The grid (dotted lines) has a box side 
length of 0.77 mm. The dotted and dashed lines are shown to 
guide the eye. (b) In the side-view cross-section (defined by 
\x — (x) | < 0.3 mm), the vertical strings of particles are clearly 
seen. The shadowed circle in (a) and shadowed rectangle in (b) 
indicate the position of an identified two-particle string. The 
trajectories of these particles during 10 s are shown in (c). The 
solid circle and solid line (open circle and dashed line) indicate 
the positions of the upper (lower) particle in the doublet. The 
doublet rotates clockwise as sketched in the inset, (d) Rela- 
tive horizontal distance between the doublet components as a 
function of time. The dotted line indicates the tracking limit 
5r < 0.2 mm. The 64-particle cluster was driven by a clockwise 
rotating electric field with a frequency of 5 kHz. 



Two or more particles can be assumed to form a string 
if they are fairly close to each other for a considerable 
amount of time that is longer than the "fly- by" time|^] 

The left panels in Fig. [3] show the histograms of the 
lifetimes of various string compositions observed in our 
experiments. The single particle lifetimes are also shown 
to emphasize the role of the interstring association and 
dissociation processes. The right panels in Fig. [3] show 
the histograms of the number of particles per string. It 
can be seen that doublet string is the most probable state 
for string configuration, and the most long-living as well. 

The observed string structure of 3D clusters agrees with 
the dust dynamics simulation of Ref. HI] . In the simula- 
tion, the clusters had a structure of nested spherical shells 
at M = 0, a state where vertically aligned pairs of parti- 



7 The fly-by time is oc A/v, where A is the interparticle separation 
and v is a characteristic particle speed. 
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Fig. 3: (left panels) Histograms of the string lifetimes At shown 
separately for strings composed of 2 to 5 particles, as well as 
single particles, (right panels) Histograms of the number of 
particles per string. Doublets and triplets are the most prob- 
able compositions of strings in a rotating cluster. The experi- 
ment was performed with a 64-particle cluster driven by a ro- 
tating electric field with a frequency of 5 kHz (rotating cluster, 
top panels) and 1 kHz (stationary cluster, bottom panels). 



cles form and move through a liquid complex plasma at 
M > 0.16, a similar "triplet state" at M > 0.33, and 
finally well aligned strings of particles at M = 1. Our ex- 
periment with estimated M = 0.6 — 0.7 fits well the above 
classification. 

String lattices. Driven clusters in our experiments 
possessed large-scale symmetry. The observed cylindrical 
shell structure had either circular symmetry [indicated by 
dashed circles in Fig. |2ja)] or tetragonal lattice symmetry 
[indicated by dotted lines] . Both states are metastable and 
the actual symmetry of the cluster changed from one to 
the other while the cluster rotated - a behavior also found 
in colloids KL8] and in quantum states in superconductors 

EH- 

String length. The spatial distribution of strings in 
a DDC is not uniform, see Fig. ga). Two main effects 
can be seen immediately: (i) the discreteness of the string 
length distribution is well pronounced and (ii) surprisingly, 
longer strings are located predominantly not in the mid- 
dle of the cluster. Another surprising fact is that on av- 
erage the triplets are relatively more compact than dou- 
blets, see Fig. gb). The overall length of a string com- 
posed, e.g. of three particles may be longer than that 
composed of only two particles. However, on average the 
three-particle strings are internally more compact, and 
this is a rather noticeable effect; the string shrinkage is 
~ 30%. (Note that 2D plasma crystals are nearly in- 
compressible in the elasticity sense [34].) The average 
particle separation for doublets is 0.87 ± 0.10 mm, while 
for triplets it is 0.78 ± 0.10 mm. For longer strings the 
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Fig. 4: String length and in-string interparticle distance, (a) 
Radial distribution of the string length. The cluster center of 
mass is located at the origin, (b) Histogram of the in-string 
interparticle distance. The string composition is color coded 
as follows: (black) all events, (red) doublets, (green) triplets, 
(yellow) quartets, (blue) quintets. The experiment was per- 
formed with a 64-particle clockwise rotating cluster driven by 
a rotating electric field with a frequency of 5 kHz. 

statistics is poor. Taking into account all events shown 
in Fig. [4^b), a mean interparticle distance in the vertical 
direction A vert = 0.83 ± 0.10 mm can be estimated. 

As a whole, the cluster consists of an outer spheroidal 
shell enclosing a cylindrically-symmetric network of mu- 
tually interacting strings. The globally spheroidal shape 
of the DDC's is most probably determined by the external 
confinement. 

Particle dragging. — We observed that in a particle 
pair in a rotating cluster, the lower particle lags behind 
the upper one, or is dragged by it, as shown in Fig.[2|c). A 
convenient way of studying this effect is to plot the distri- 
bution of azimuthal angles Acp of particle pairs, see Fig. [5] 
(bottom panel). In the stationary cluster (black line), A(f 
has a broad distribution roughly centered at zero. In the 
counterclockwise rotating cluster (red line), an additional 
peak develops at ~ —tt/2. In the clockwise rotating clus- 
ter (blue line), the peak shifts to w 7r/4. These peaks at 
Aip 7^ indicate the particle dragging effect. The reason 
why the peak in the latter case appears at « 7r/4 instead 
of the expected value of tt/2 is not clear. 

Observed particle dragging is an important effect, which 
directly demonstrates the presence of in-string attractior^\ 
In our experiments, the driving torque on the cluster de- 
pends on the vertical coordinate [22]. A dragging is then 
to generate an effective coupling mechanism in response 
to the driving torque inhomogeneity. In other words, the 
transverse component of the attraction force to the ion fo- 
cus F w is to compensate for the moment deficit AM.. F w 
evidently could not exceed the friction force Ff r , that is: 



j- id 



AM 



< F h = mjefiU) c R, 



(2) 



where R is the distance to the cluster center, uj c is the clus- 
ter rotation speed. The relationship (J2| yields an upper 

8 To some extent, this resembles the laser-dragging experiment of 
Ref. fi2l. 
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Fig. 5: Probability densities of (top) polar and (bottom) az- 
imuthal A(p angles for particle strings in a 3D cluster. The 
probability densities are calculated from the respective distri- 
butions by normalizing the area under the curve. They are 
shown for three situations: counterclockwise rotation (red line, 
rotating electric field at 5 kHz), clockwise rotation (blue line, 
5 kHz), and no rotation (black line, 1 kHz). The angle was 
measured between the vertical z-axis and the (imaginary) line 
connecting the upper and lower particles of a string, (top in- 
set) P(0) is fitted well with Eq. [I] where A = 188.2 ± 3.7, 
(60 2 ) = 0.0104 ± 0.0002. (bottom inset) Aip is defined as the 
relative angular position on the (x,y) plane of the upper (U) 
and lower (L) particles of a string. C is the cluster center of 
mass. The range of Aip was extended to [— 2tv, 2ty] to take into 
account two turns made by a particle string. 



estimate F w < 10 fN for our parameter^] 

It is natural to introduce a "spring constant" k w of this 
wake-mediated interaction: 



w 

W)' 



(3) 



where (Sr) is the mean relative radial displacement of the 
dragged particle. For instance, for the doublet traced in 
a rotating cluster shown in Fig. [2jc),(d) (Sr) = 0.07 ± 
0.04 mm, which yields k w « 900 eV/mm 2 . 

Vertical alignment. — The strings that form inside 
DDC's in our experiments are highly aligned vertically. 
This is clearly seen in the strings' polar angle distribution, 
Fig. [5] (top panel) and is also corroborated by large values 
(about 80 — 90%) of the vertical order parameter 35 . This 
is consistent with the results of dust dynamics simulations 



of Ref. [14] for our estimated ion flow Mach number of 
M = 0.6 -0.7. 

The polar angle probability densities happened to be 
similar in shape regardless of whether the DDC is clock- 
wise or counterclockwise rotating, or even static. Such 
kind of universality is actually not surprising because the 
rotational degrees of freedom, in particular the polar angle 
distribution, must obey thermal statistics: 

P{6) = AO exp (- jj^j , 59 = 9- (9) , (4) 

where A is a coefficient and (SO 2 ) is the mean squared 
variation. The value (SO 2 ) ~ 0.01 gives a fairly good fit 
for all cases shown in Fig. [5] By virtue of the fluctuation- 
dissipation theorem, it is straightforward to obtain the 
kinetic temperature T of the particle chaotic vibrations in 
the DDCE3 

U *° (SO 2 ) ^ 3.2 eV. (5) 



-h A 2 



While being approximate, the relationship ([5| provides 
a fairly good estimate of T = (2/3) (rm; 2 /2), where the 
mean kinetic energy of particles (mv 2 /2) « 4 eV was ob- 
tained directly by particle tracking technique. 

String-string interaction. — The dynamics of 
paired particles in strings (as well as of isolated pairs) 
can be successfully used to study the in-string and string- 
string coupling constants. Evidently there are different 
time scales involved in this analysis. The short time scale 
evolution At < 0.1 s can hardly be addressed properly in 
view of the strong pixel locking effect (for details, see e.g. 
Refs. 36,37 ). The longest time scale of the order of a few 
rotation periods At > 20 s is difficult to study since one 
needs to operate with enormously large amount of data. 
(In the given experiment we limited ourselves to At = 10 s, 
which is approximately two thirds of the rotation period.) 
The remaining intermediate time scale At ~ 1 s is the 
most important for our purposes. It allows one to study 
the evolution of the global structure of the cluster, which 
is determined predominantly by the mutually interacting 
strings. 

The interaction between strings includes direct particle 
exchange [9] and electrostatic interaction. From the spec- 
trum of string deformations, the spring constant of the 
string-string interaction can be estimated. As a standard 
Fourier analysis revealed, the time scale of the string os- 
cillations corresponds to frequencies of 0.5 — 3 Hz, which 
are easy to study. Such a low frequency domain of string 
deformation can be successfully described assuming a dy- 
namical balance between the friction force and the local 
force deforming the string. Assuming that the deforma- 
tion is in the elasticity domain, with k* being the Hooke's 
spring constant, one can conclude that 



27T7eff/def^, 



(6) 



9 To compare, the force scale of Yukawa-interacting particles in a 
sting is about 200 — 300 fN, the energy scale is about 300 — 500 eV. 



10 Here, only as a kinetic average, since DDC's are not in equilib- 



p-5 



where /def is a characteristic frequency of the low- 
frequency deformation. For /def ~ 1 — 2 Hz the relation 
^ yields k* ~ 600 - 1200 eV/mm 2 . It is noteworthy that 

To conclude, 3D complex-plasma clusters driven exter- 
nally by a rotating electric field were experimentally ob- 
served to be highly structured. They showed ion-flow 
aligned string lattices of a cylindrical symmetry. A fully 
3D holographic particle tracking diagnostic allowed us to 
thoroughly explore the statistical as well as dynamical pa- 
rameters of the clusters' string structure. Rotating electric 
field was shown to be a useful manipulation tool for 3D 
complex-plasma clusters. 

* * * 
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